Chaotic instantons and enhancement of tunneling in double-well system with 
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Chaotic instanton approach allows to describe analytically the influence of the polychromatic 
perturbation on quantum properties of nonlinear systems. Double well system with single, multiple 
and polychromatic kicked perturbation is considered in the paper to compare quantitative analytical 
predictions with the results of numerical calculations. Chaotic instantons are responsible for the 
enhancement of tunneling far away from the exact (avoided) level crossings in framework of the 
approach used. Time-independent averaged system is used for regular approximation of the chaotic 
instanton solution in order to take into account its contribution to the ground quasienergy doublet 
splitting. Formula for the ground quasienergy levels' splitting was derived averaging trajectory 
action in the stochastic layer in the framework of chaotic instanton approach. Results of quasienergy 
spectrum numerical calculations and simulations for tunneling dynamics are in good agreement with 
the obtained analytical predictions. 

PACS numbers: 03.65.Xp, 03.75.Lm, 05.45.Mt 



I. INTRODUCTION 

Investigation of the influence of small perturbation 
on the behavior of the nonlinear dynamical systems 
continues to attract permanent interest for several last 
decades [H-[3|- The connection between the semiclassi- 
cal properties of perturbed nonlinear systems and purely 
quantum processes such as tunneling is a reach rapidly 
developing field of research nowadays [!, Q . Our insight 
in some novel phenomena in this field was extended dur- 
ing the last decades. The most intriguing among them 
are the chaos assisted tunneling (CAT) and the closely 
related coherent destruction of tunneling (CDT). 

In particular, the former is an enhancement of tunnel- 
ing in the perturbed low-dimensional systems at small 
external field strengths and driving frequencies 
This phenomenon takes place when levels of the regu- 
lar doublet undergo an avoided crossing with the chaotic 
state pil H3- At the semiclassical level of description one 
considers tunneling between KAM-tori embedded into 
the " chaotic sea" . The region of chaotic motion affects 
tunneling rate because, compared to direct tunneling be- 
tween tori, it is easier for the system to penetrate pri- 
marily into the chaotic region, to travel then along some 
classically allowed path and finally to tunnel onto another 
KAM-torus [H,[l3|. The latter, CDT phenomenon, is a 
suppression of tunneling which occurs due to the exact 
crossing of two states with different symmetries from the 
tunneling doublet [lij]. In this case the tunneling time 
diverges which means the total localization of quantum 
state on the initial torus. 

CAT phenomenon as well as CDT were experimen- 
tally observed in a number of real physical systems. The 
CAT observation between whispering gallery-type modes 
of microwave cavity having the form of the annular bil- 
liard was reported in Ref. [Hj]. The same phenomenon 
for ultracold atoms was experimentally investigated in 



Refs. [H, The study of the dielectric microcavities 
provided evidences for CAT as well [l8[ . Both CAT and 
CDT phenomena were observed in two coupled optical 
waveguides [H, ■ Recently experimental evidence of 
coherent control of single particle tunneling in strongly 
driven double well potential was reported in Ref. [21]. 

The most common methods which are used to investi- 
gate the interplay between semiclassical properties of per- 
turbed nonlinear systems and quantum processes are nu- 
merical methods based on Floquet theory [!, [lj|, [22| and 
Random Matrix Theory (23[. Among other approaches 
we would like to mention the scattering approach for bil- 
liard systems [24], HE} and approach based upon the pres- 
ence of a conspicuous tree structure hidden in a compli- 
cated set of tunneling branches [26H28l | . 

In this paper we will consider the original analytical 
approach based on instanton technique. Enhancement of 
tunneling in system with external force in framework of 
this approach occurs due chaotic instantons which ap- 
pear in perturbed case. This approach was proposed 
in Refs. [29l - l3l| and used in Ref. 32]. Chaotic in- 
stanton approach will be developed further here using 
averaged time-independent Hamiltonian and exploited 
for description of the enhancement of tunneling in the 
polychromatically kicked double well system. Previously 
polychromatic perturbation was investigated numerically 
only [i^]. The main purpose of the present study is to 
prove the ability of developed chaotic instanton approach 
to give quantitative analytical description of tunneling in 
polychromatically perturbed systems well agreed with in- 
dependent numerical calculations based on Floquet the- 
ory. It will give additional support and pulse for the 
further development of analytical methods to investigate 
tunneling phenomenon in quantum systems with mixed 
classical dynamics. Alternative approach based on quan- 
tum instantons which are defined using an introduced 
notion of quantum action was suggested in Refs. [33l l34j] . 
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Another instanton approach was developed recently for 
description of ordinary, dynamical and resonant tunnel- 
ing in various nonperturbed systems [35]. Analytical ap- 
proach to describe tunneling in perturbed systems based 
on nonlinear resonances consideration was developed in 
Refs. [3^, [13] ■ A theory for dynamical tunneling process 
using fictitious integrable system was recently given in 
Refs. [HSi]. 

Double well potential is a model which is convenient 
to use for tunneling analysis. This system is well stud- 
ied in the nonperturbed case, e.g. on the base of in- 
stanton technique (Hill or WKB method Double 
well potential is often used for description of processes 
which occurred in wide range of real physical systems: 
such as flipping of the ammonia molecule (43| . transfer 
of protons along hydrogen bonds in benzoic-acid dimers 
at low temperatures (4J, 1451 ] and macroscopic quantum 
coherence phenomena in superconducting quantum inter- 
ference devices and nanomagnets [49|, l50| ■ Per- 
turbation in this paper is regarded in the form of kicks. 
One of the attractive features of this type of perturba- 
tion is the extensively-investigated simple quantum map 
which stroboscobically evolves the system from kick n to 
kick n + 1. Kicked systems are recently used for experi- 
mental realization of a such novel concept as a quantum 
ratchet [5ll [52| . Double kickedperturbation was investi- 
gated experimentally in Ref. [531 ] . 

The paper is divided into several sections. Chaotic in- 
stantons are analyzed using averaged time-independent 
Hamiltonian of the kicked system in section [TTJ Re- 
sults obtained by means of the averaged Hamiltonian 
are used in section IIIII to derive analytical formula for 
lowest quasienergy doublet splitting dependence on per- 
turbation parameter in single kicked system. Numerical 
calculations are performed to check the validity of this 
formula in the section IIVI Multiple and polychromatic 
kick perturbations are considered in the sections [V] and 
[VTl respectively. 



II. CHAOTIC INSTANTONS IN KICKED 
DOUBLE- WELL POTENTIAL 

Hamiltonian of the particle in the double- well potential 
can be written down in the following form: 

2 

Hq = - h a Q x 4 - a 2 x 2 , (1) 

Zm 

where m - mass of the particle, ao, a 2 - parameters of the 
potential. We consider the perturbation of the kick-type 
and choose it as follows: 

+00 

V per = eTx 2 Ht-nT), (2) 



coordinate was chosen in the form of x 2 in order to pre- 
serve spatial symmetry in the perturbed system. Hamil- 
tonian of this system is 

H = H + V per . (3) 

Now we implement Wick rotation (t — > —it) and define 
Euclidean Hamiltonian 

+00 

U E = U E - eTx 2 S(r-nT), (4) 

Tl — — OO 

where T~L E - nonperturbed Euclidean Hamiltonian which 
is given by 

r? 

uE V 4,2 /r\ 

H = a x + a 2 x . (5) 

Zm 

Euclidean equations of motion of the particle in the 
nonperturbed double-well potential (e = 0) have a well 
known solution - instanton. This solution is used for 
calculation of the ground energy splitting in the system 
without perturbation [!(], 5l| and explains the rate of 
the tunneling process in it. Another solutions of the Eu- 
clidean equations of motion besides ordinary instanton 
are required to explain dynamical tunneling in perturbed 
system. Perturbation destroys the separatrix and some 
trajectories in its vicinity go to infinity. Narrow stochas- 
tic layer is formed nearby the nonperturbed separatrix 
due to the perturbation. "Chaotic instanton" appears 
in this layer. Chaotic instanton is the closest to the de- 
stroyed nonperturbed separatrix trapped trajectory (see 
figure [1} . Thus it plays a dominant role in tunneling in 
perturbed system. 

Now let us construct the averaged time-independent 
Hamiltonian for the double well system with the pertur- 
bation of the kick-type. It is calculated using Euclidean 
Hamiltonian ((4]) 

H*=±;J H E dt = ^-a x A + (a 2 -e)x 2 . (6) 

Hamiltonian ([6]) is coincided with nonperturbed 
Hamiltonian ([3]) when parameter a 2 is replaced by a 2 . 
The latter is defined as follows: 

a 2 (e) = a 2 - e. (7) 

In contrast to the kicked system ((4]) averaged Hamilto- 
nian ^ is time-independent. Since the form of the po- 
tential is changing, there should be restriction for the 
perturbation strength variation. This restriction follows 
from our assumption that ordinary instanton approach 
should be valid for the averaged potential. Condition 
for ordinary instanton approach applicability [4ll | can be 
written down for instanton action in averaged potential 



where e and T are perturbation strength and period, re- 
spectively, t - time. Dependence of the perturbation on 
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FIG. 1: Phase space of the system with averaged Hamilto- 
nian with parameter 52 (e = 0.02). Separatrix in this system 
(thick solid line) and in the nonperturbed system (dashed 
line) are shown in the figure. Comparison of the particle clas- 
sical motion on one period of the perturbation in averaged 
(thick solid lines) and kicked (thin solid lines) systems from 
the set of initial conditions (thick points) are shown in the 
inset (a). Inset (b) shows a stroboscopic plot of chaotic tra- 
jectory of the kicked particle near the turning point of the 
separatrix in averaged system. 



where ujq - oscillation frequency near the bottom of the 
wells, Si ns t = 2^/m a^ 2 / (3 ao) - nonperturbed instanton 
action. 

Ground doublet quasienergy splitting (A77) in the 
kicked system in the framework of our approach is 
expressed in terms of chaotic instanton action (S c h) 
through the formula which is similarly to : 



Ar/ = 2u>ov - 




) • <'->• 



72 S c h 



where chaotic instanton action can be calculated by av- 
eraging the nonperturbed trajectory action (S(E)) over 
energy for stochastic layer width 



S c h — 



1 



AH 1 . 



S{E)dE = 



1 



AU{ 



where we have made a transformation to the integral over 
the energy difference £ = Ei ns t—E in last expression. Us- 
ing nonperturbed trajectory action expansion near the 

separatrix S(E) = irJ(E inst - £) S mst - a £ ex- 
pression for chaotic instanton action can be calculated 
directly. Here a — (1 + 181n2)/6 is a numerical coeffi- 
cient. Thus, for chaotic instanton action we obtain 



Thus, for restriction considered wc obtain 



£ $ tmax — 0,2 



3^3 ajf. 



(9) 



Figure Q] shows a possibility to describe properties of 
the classical motion in the kicked double well system in 
Euclidean time using averaged system ©. This Hamil- 
tonian was used to analyze the perturbed system phase 
space in Euclidean time and to construct an approxima- 
tion for chaotic instanton solution. This approximation is 
a separatrix in the averaged model. Using it's properties 
we obtain the formula for stochastic layer width 



2a 



(10) 



Expression (|10p will be used in the following section 
in order to obtain analytical formula for the lowest 
quasienergy doublet splitting dependence on the pertur- 
bation strength in the kicked system. 



III. GROUND DOUBLET QUASIENERGY 
SPLITTING FORMULA 



The lowest doublet energy splitting in two loop ap- 
proximation in the nonperturbed double well potential is 
the following (see |54{ and review |41j): 



Sch — So — —t — A%^ 

4 V °2 



(13) 



Now we can write down analytical formula for 
the ground quasienergy levels splitting using expres- 
sions (HFJ, (HIJ), (dU) and fl3|: 



Arj{e) = AE e* 



where 



k = 



a y/m 0,2 
4 a 



(14) 



(15) 



Tunneling period in the kicked double well potential is 
expressed in terms of ground quasienergy levels splitting 



2tt 
At]' 



(16) 



Increasing of the perturbation parameter gives exponen- 
tial rise to ground quasienergy splitting and to the tun- 
neling frequency (ft un (e) = A 77(e)). The last exponen- 
tial factor in the expression (ITU is responsible for the 
tunneling enhancement in the perturbed system. In non- 
perturbed case formula (| 14[) coincides with the expres- 
sion (fTTj) . Formulas (|14p and (|16p will be checked in 
numerical calculations in the next section. 



IV. NUMERICAL CALCULATIONS 



A£o = 2 cj 



st exp 



Si', 



71 1 

72 Sint 



, (11) 



For the computational purposes it is convenient to 
choose the eigenvectors of harmonic oscillator as the basis 
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vectors. In this representation matrices of the Hamilto- 
nian ([1} and the perturbation are real and symmetric. 
They have the following forms (n > m): 



h,c [n+^j + ^(^ga (2m 2 + 2m + 1) 



ttO _ r 

- a' 2 (2m + l))] 



+ S m +2n X (ffao(2m + 3) - a 2 ) y/(m + l)(m + 2) 



S m +A n W + + 2)(m + 3)(m + 4), 



where <? = h/mu> and a' 2 = 02 + mcj 2 /2, is Planck con- 
stant which we put equal to 1, uj - frequency of the basis 
harmonic oscillator which is arbitrary, and so may be ad- 
justed to optimize the computation. We use the value 
u> = 0.2 with parameters m = 1, ao — 1/128, a-i = 1/4 
in most of calculations which are chosen in such a way 
that nonperturbed instanton action is large enough for 
energy splitting formula for nonperturbed system to be 
valid and not too big in order to decrease errors of numer- 
ical calculations. The matrix size is chosen to be equal 
to 200 x 200. Calculations with larger matrices give the 
same results. System of computer algebra Mathcmatica 
was used for numerical calculations. 
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FIG. 2: Quasienergy spectrum for the ten lowest average en- 
ergy levels. Thick lines - doublet with the minimal average 
energy. 



We calculate eigenvalues of the one-period evolu- 
tion operator ( e ~ iH oT/2 e -ieTx e -iH T/2^ an ^ obtain 
quasienergy levels (rjk) which are related with the evo- 
lution operator eigenvalues (A&) through the expression 
rjk = i In Afe /T. Then we get ten levels with the lowest 
one-period average energy. We obtain these levels us- 
ing the formula (vi\H av \vi). Here H av is the averaged 
Hamiltonian in Minkowski space, \vi) are the eigenvec- 
tors of the one-period evolution operator. The depen- 
dence of ten lowest levels' quasienergies on the strength 
of the perturbation for the model parameters mentioned 



above is shown in the figure [2] Quasienergies of two lev- 
els with the minimal average energy are shown by thick 
lines. They are too close to each other to be resolved in 
the figure due to very small splitting. 
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FIG. 3: Quasienergy splitting as a function of the strength 
of the perturbation for different values of model parameters 
ao (a) and 02 (b). Lines - analytical formula (|14l) . points - 
numerical results. 

Performed numerical calculations give the dependence 
of the ground quasienergy splitting on the strength of the 
perturbation for different values of model parameters ao 
(figEHa)) and 112 (ngj3jb)). We fix parameter a 2 = 1/4 
for figure [3Ja) and ao = 1/128 for figure [3fb). Re- 
sults of numerical calculations are plotted in the figure [3] 
by points. Axis Arj is shown in logarithmic scale. Ob- 
tained dependencies are exponential as it was predicted 
by chaotic instanton approach and obtained analytical 
formula (fT4"j) . 

Analytical results are plotted in the fig- 
ures |3] (a) and (b) by straight solid lines. Numerical 
points lie close to these lines. The agreement between 
numerical calculations and analytical expression is good 
in the parametric region considered. 

Now lets perform numerical simulations for the tunnel- 
ing process in the kicked double well system and check an 
applicability of the formulas (fT4"|) and (fTBT) for this pro- 
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confirms analytical assumptions mentioned above. 





FIG. 4: Quantum mechanical tunneling in kicked double 
well potential. Perturbation parameters: (a) e = 0, (b) 
T = 2tt/4, e= 1.9 ■ 10~ 2 . 



cess. For this purpose we regard the double well poten- 
tial (H} with parameters m — 1, a — 1/128, a 2 = 1/4 
and the same basis vectors as for previous calculations. 
We take a symmetric superposition of two lowest nonper- 
turbed states as a initial wave packet. These packet is 
localized in left well of potential. Numerical simulations 
we provide by multiplying initial wave function by one 
period evolution operator. The results of numerical sim- 
ulations for the two values of the perturbation strength 
are shown in the figure |U The dependence of the local- 
ization probability of the wave packet on the coordinate 
and time is presented in figures. Minima of the nonper- 
turbed double well potential ((J) are situated in points 
x = —4 and x = 4. Tunneling between these points in 
nonperturbed system is demonstrated in the figure 21a) . 
Evolution of the initial wave packet in perturbed case is 
shown in the figure HJb). Perturbation parameters for 
these simulations are T = 2-7r/4 and e = 1.9 ■ 10 -2 . They 
are chosen in such a way to speed up a tunneling in two 
times in comparison with nonperturbed system. Figures 
Ufa) and BJb) demonstrate this enhancement. Fourier 
analysis of the dependence of the localization probability 
of the wave packet in left well on time in perturbed case 
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FIG. 5: Comparison error of the analytical formula ()14|) with 
results of numerical calculations in percents. 

In order to check applicability of the developed ap- 
proach we carry out a series of the numerical calculations 
for wide range of the perturbation parameters. Results of 
the analysis is performed in figure [5] Region of quantita- 
tive agreement between analytical and numerical results 
is shown by black color in the figure. There are two re- 
strictions of the developed approach applicability. The 
first one is that model parameters should be far away 
from the exact (avoided) level crossings. Thus we have 
the restriction for the perturbation period (T <C 27r/wn, 
where ujq - oscillation frequency near the bottom of the 
wells). Another restriction for analytical predictions is 
a condition for ordinary instanton approach applicability 
which imply the maximum for perturbation strength (0) . 
This maximum is shown in the figure by dashed line. The 
two restrictions mentioned above explain accurately the 
figure [5l 

Inverse sign in the expression of the perturbation @ 
will induce exponential suppression of tunneling in the 
system. It can be demonstrated numerically as well. 

V. MULTIPLE KICK PERTURBATION 

Lets consider double-well system with multiple kick 
perturbation. Hamiltonian of this system is the follow- 
ing: 

2 +oo 

H = — + a .x 4 - a 2 x 2 + e{Tx 2 V S(t - nT) 

n= — oo 

+ OC 

+ e 2 Tx 2 S(t + AT-nT), (17) 

n— — oo 

where t\ and e 2 are strength of two perturbations, T is 
the period for both kicking sequences, AT - shift between 
these sequences. 
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Averaged time-independent Hamiltonian for the sys- 
tem under investigation is given by 



1 

H a v — y 



Hdt 



P , 4 

2m 



{a 2 -e 1 -e 2 )x 2 . (18) 



Using the last expression we can rewrite restriction for 
perturbation strength ([8]) in multiple kick case 



ei + £2 < = a 2 - 3 ^3 a 2 



2/3 



(19) 
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where coefficient fc is defined using expression (I15I) . 

Obtained analytical formula (|2T)]) is checked in numer- 
ical calculations (see figures [6] and [7]) . One period evolu- 
tion in case considered has the following form: 

(21) 

Figure [6] shows that formula (|20l) can be used for the 
description of the ground quasienergy splitting depen- 
dence. The applicability parametric region of the for- 
mula is demonstrated on the figure [7] Restriction (fT9|) is 
shown by dashed line. 



VI. POLYCHROMATIC PERTURBATION 

Finally we will consider double kick system with dif- 
ferent values of the perturbation period. System Hamil- 
tonian is given by 



tt P , 4 2 , rr 

H = - h a^x ~ a 2 x + eilix 

2m 



^T 2 x 2 S(t-nT 2 ) 

n= — oo 



3 £ 

n— — oo 



S(t - nTi) 



(22) 



FIG. 6: Quasienergy splitting as a function of the strength of 
the second perturbation for different values of the first one. 
Perturbation period T = 1, kicks shift dT = 0.4. Lines - 
analytical formula, points - numerical results. 



where T\ and T 2 are periods of two perturbations. 

Averaged Hamiltonian can be calculated by averaging 
of the perturbed Hamiltonian over time for less common 
multiple Ti cm of two periods T\ and T 2 



Chaotic instanton approach applicability error in percents 
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FIG. 7: Comparison error of the analytical formula (|20|) with 
results of numerical calculations in percents. Perturbation 
period T — 1, perturbation strength ei = 0.01. 

Using expressions for quasienergy splitting (TT2"]) . 
chaotic instanton action (TT3")) and averaged Hamilto- 
nian (ITS)) analytical formula for ground quasienergy split- 
ting in the multiple kicked double well potential is ex- 
pressed in terms of perturbation strength values 



Arj(e) = AE e kitl+e2 \ 



i 



Ti 



Icm JO 



Hdt = — — h anm 4 — (02 — e\ — e 2 ) x 2 . 
2m 

Obtained averaged Hamiltonian is the same as for 
monochromatic multiple kicked system (|18p . Thus ana- 
lytical formula for ground quasienergy splitting will have 
the form (l20l). 



Chaotic instanton approach applicability error in percents 
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FIG. 8: Comparison error of the analytical formula (|20l) with 
results of numerical calculations in percents. Perturbation 
periods Ti = 0.6 and T 2 = 1. 
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Chaotic instanton approach applicability error in percents 
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FIG. 9: Comparison error of the analytical formula (|20p with 
results of numerical calculations in percents. Parameters of 
the first perturbation T\ — 1 and ei = 0.01. 



To calculate quasienergy levels we construct evolution 
operator U (TJ cm ) for period of time T\ cm in analogy with 
expression (j!?Tj) . Subsequent steps of numerical calcu- 
lations are identical to the algorithm presented in the 
section [IV] 

Maps of the approach applicability for aperiodic kicked 
system are presented on figures [5] and |H1 There is a good 
agreement between analytical and numerical results in 
a wide range of perturbation parameters. Chaotic in- 
stanton approach applicability restrictions T <C 2tt /ujq 
and (fT5|) explain the obtained figures. The last restric- 
tion is shown by dashed line on figures. 

VII. CONCLUSIONS 

Chaotic instanton approach allows to describe analyt- 
ically the influence of the polychromatic perturbation on 



quantum properties of nonlinear systems. Double well 
system with single, multiple and polychromatic kicked 
perturbation is regarded in the paper to compare quanti- 
tative analytical predictions with the results of numerical 
calculations. 



Chaotic instanton is the solution of the Euclidean 
equations of motion of the perturbed system. This con- 
figuration is responsible for the enhancement of tunneling 
far away from the exact (avoided) level crossings. Time- 
independent averaged system is used for regular approx- 
imation of the chaotic instanton solution in order to take 
into account its contribution to the ground quasienergy 
doublet splitting. Formula for the ground quasienergy 
levels splitting was derived averaging trajectory action 
in stochastic layer in the framework of chaotic instanton 
approach. This formula predicts exponential dependence 
of the ground doublet splitting on value of the perturba- 
tion strength. 



Numerical calculations for quasienergy levels depen- 
dence on value or values of single, multiple and polychro- 
matic perturbation strength and simulations for tunnel- 
ing dynamics are performed to check the validity of the 
obtained analytical formulas. Results of numerical cal- 
culations for the quasienergy spectrum confirm the ex- 
ponential dependence of the ground splitting on value of 
the perturbation strength for single perturbation or sum 
of values in multiple kicked case. They are in good agree- 
ment with the derived analytical formulas (fL4")l and (|20l) . 
Simulations of the tunneling dynamics in the kicked dou- 
ble well system demonstrate exponential tunneling en- 
hancement as well. Applicability of chaotic instanton 
approach was tested in a series of numerical calculations. 
Sufficiently wide range of perturbation parameters was 
found suitable for developed approach application. 
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